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We analyze the consequences of a first-order (in the effective Kondo coupling) inter-impurity 
correlation present in the Two-Impurity Anderson Hamiltonian. This correlation generates a pair 
of Kondo Doublet states which are responsible for the first (higher energy) Kondo stage of these 
systems. Applied to systems of quantum dots in a semiconductor heterolayer, these Doublets predict 
a mainly ferromagnetic dot-dot spin correlation which is different than the expected from the second- 
order RKKY interaction. In order to effectively screen both dots, a parallel alignment of the dot- 
spins is favored by the interaction. 
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The Anderson Impurity Hamiltonian describes the in- 
teraction of a localized electron state (with a high internal 
Coulomb interaction) with a set of uncorrelated-extended 
orbitals. Primary used to study the behavior of magnetic 
impurities in a metallic host[l|, it also applies to 

systems of quantum dots (QD) connected to/through a 
low-dimensional electron gas (LDEG), as has been cer- 
tified by the observation of the Kondo effect in a single 
quantum dot, a pioneer work of Goldhaber-Gordon and 
coworkers^. QD states with a well-defined number of 
electrons tend to be stabilized by the Coulomb blockade 
effect. When an odd number of electrons is stable in the 
dot and its total spin is 1/2, the hybridization coupling 
with the LDEG gives rise to the usual Kondo effect. 

A logical step in the development of nanoscopic elec- 
tronics is to considerer systems with several QD0, 0, 
0, 0, in fact arranges of QD has been proposed with 
the localized QD-spin (Si) state playing the role of a 
Quantum- bit. For such devices to be useful, a detailed 
knowledge of the behavior of the QD-spin correlations 
(Sj.Sj) is needed. This correlation is expected to be con- 
tained in the two-impurity Anderson model. 

The two-impurity Anderson Hamiltonian has been 
the subject of many theoretical studies, ranging from 
perturbation theory and narrow band approximation 
l[ p]. [To [ [ill to numerical renormalization group analysis 
HE13E1II1II3 and conformal field theory [3. An 
other theoretical approach to Anderson Impurity systems 
has been the use of variational wave functions (VWF) 
H 13 El Ell In Ref. H3 a set of VWF equations that 
describes the lower energy state of the system has been 
numerically anal yze d, these equations were analytically 
solved in Ref . [H |23|] . These VWF are particularly useful 
to analyze internal correlations between the components 
of the system because they are constructed following such 
correlations. 

In what follows we analyze the consequences of a very 
simple inter-dot correlation that generates a pair of corre- 
lated Kondo Doublets states. Surprisingly, these Kondo 



Doublets generates a mainly ferromagnetic QD-spin cor- 
relation without resort to the second-order RKKY inter- 
action. 

The Anderson Hamiltonian for magnetic impurities (or 
QD) connected to a metallic electron gas (a LDEG) is the 
sum of the band Hb(ek) and effective QD Hd(—Ed,U) 
Hamiltonians plus the hybridization term Hy(V) 

H = J2 e k c\ a c ka + v ^(e 1 k - r J d] a c ka + h.c.) 

ha jka 
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where the fermion operator Cka(dj a ) act on the conduc- 
tion band fc-state (on the QD state situated at rj). Single 
state energies e k ,—Ed are refereed to the Fermi energy 
(Ep), v = V/y/N is the d — c hybridization V divided 
by square root of the number of cells in the metal (N) . 
In the Kondo limit, the case analyzed in this paper, the 
QD level is well below the Fermi energy (— Ed <C 0), 
and they can not be doubly occupied due to the strong 
Coulomb repulsion in them (U — > 00). We renormalize 
the vacuum (denoted by \F)) to be the conduction band 
filled up to the Fermi energy and we make an electron- 
hole transformation for band states below the Fermi level: 
61 = c-k.-o for \k\ < kp; with the above renormaliza- 
tion the hole excitation energies are explicitly positive. 

We consider here the two QD case, one QD at — r/2 
and the other at r/2, over the a;- axis. We use a "ket" 
notation for the QD population, the first symbol indi- 
cates the status of the left QD and the second one the 
status of the QD on the right, e.g. |J.T) = dLd^JF), 

|0T)s4 T |F), |00) = \F), etc.. 

The two relevant Hamiltonian parameters in the 
regime detailed above result to be the effective Kondo 
coupling J n = n V 2 / Ed (n being the density of band 
states at the Fermi level) and the inter-Dot distance r. 
The half band-width D sets the energy scale. 

Correlations between the QD are the consequence of 
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closed electron-paths that involve both QD, paths driven 
by Hy in configurational space. The four lowest energy 
states are the ferromagnetic triplet (FM, Iff),...; for 
which (Sl.Sr) — 1/4 ) and the antiferromagnetic singlet 
(AF, Ut) - I Tl); with (S L .S R ) = -3/4), with an energy 
—2Ed- A well known second order (in J„, four in v) 
correlation is the RKKY 

process HE!! , 

which splits 

the FM triplet from the AF state. Starting from the FM 
state an electron (say the one in the left QD) hops to a 
band state k(> Uf) and then a second electron hops into 
the QD, letting a hole q(< kp) in the band. The path is 
closed by reversing these steps but at the other QD. The 
matrix element for this loop is 



<TT |i4|TT) m r = Tv 4 cos(k + q).r 



(2) 



where the plus sign holds for the AF state. This cor- 
relation path produces an energy correction =pE^(r) 
for the FM (AF) state which splits the FM and AF 
states, it is easily evaluated by perturbative0 or VWF 
methods j^. 

The proper one-electron-in-two-QD states are the sym- 
metric and antisymmetric localized spin 1/2 states 

K) = -1(4, =p 4)|F) = -L(\0a) T WO)) , (3) 

where the plus sign holds for the symmetric state \S a ), 
these states have an energy —Ed. They are the starting 
point of a first order correlation path. The hybridization 
Hy connect the \A-\) state, by promoting an electron 
from below kp to the empty QD, with the following ones 



\A Uk ) = -J bl (e- ik --/ 2 \U)+e +lk - r/2 UT» , (4a) 



V2 
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\A m ) = -J b{ L (e- ik ^ 2 |tt> +e +Jk - r/2 |TT)) , (4b) 



of energy —2Ed + e k . Closing the loops, as depicted in 
FigH the following matrix elements are obtained 

(A^\H v \A tlk ) = v 2 , (5a) 
(A^\H v \A nk ) = v 2 (l±cosk.r) , (5b) 

where the minus sign corresponds to a similar path but 
for the \S-() state. The last element depends on the inter- 
QD distance r, thus it gives a correlation between the 
status of the two QD. 

As the starting state of the loop )) has a higher en- 
ergy than the bunches of "visited" ones (l-Afo-fc)) the en- 
ergy associated with this correlation can not be obtained 
by perturbative methods 0, 0] . This situation, a single 
state connected to bunches of lower energy states, is the 
hallmark of a Kondo structure, in this case a Kondo odd 
doublet that can be analyzed by means of the following 

vwfEHH 



(6) 



|D oT ) = L4 T )+E Z(k)(\A nk ) + \A m )) . 
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FIG. 1: First order correlation paths for the \A^) QD spin 
configuration. On the left it is shown the non-correlated aa 
channel. On the right it is the interference-enhanced aa chan- 
nel, the crossed arrows correspond to the inter-QD paths. 



Proceeding as in the one impurity case|2j, the variational 
amplitude turns out to be Z(k) = l/(E Q + 2Ed — e k ), 
where E a is the energy of the doublet 

E = -E d + v 2^ - ; OF — — , (7) 

k 



E + 2E d - e k 



assuming E a = —2Ed — S (thus Z(k) = —l/(5 + e k )) 
we obtain 



E d + 5 = v 2 Y, 



2 + cos k x r 
S a + e k 



(8) 



a self-consistent equation for the inter-QD correlation 
energy S a (r). Two limits are easily found. For very 
large r the contribution from the cos term vanishes and 
thus 5 a tends to 5k = Dexp (— 1/2 J„), the one-impurity 
Kondo energy. At r = (cosfc^r = 1) one obtains 
MO) = S 3 = L>cxp(-1/3J„) » 5kW&- This is a huge 
correlation energy, for J n = 0.0723 (i.e. 5k = 0.001D) 
one obtains a tenfold increment over the one-QD Kondo 
case, given that 5% = 0.0 ID. 

For arbitrary inter-QD distances one must evalu- 
ate the sums in Eq.JSJ). One is the Kondo Integral 
I K (6) = E k (-Z(k)/N) = n a hi(l + D/5) and we call 
the other the Quantum Coherence Integral Iq(5, r) = 
S fe (- cos (k x r)Z(k)/N) which for a ID LDEG becomes 
Iq(S,r) = n (cos(y)[Ci(y) - Ci(x)] + sin (y) [Si(y) - 
Si(x)]) , where x — kpr and y = x(l + 5/D), Ci (Si) 
is the Coslntegral (Sinlntegral) function 22]. The Co- 
herence Integral has a logarithmic dependence on 5 that 
goes like Ik, note that Iq(S, 0) = Ik (5), so it is useful to 
define the coherence factor 



C Q (5,r) = I Q (5,r)/I K (6) , 



(9) 



which is a decaying oscillatory function that depends 
weakly on 5. For low dimensional or spatially confined 
LDEG, such that the main decoherence factor is the en- 
ergy width of the packet of electrons involved in the 
interaction, Cq remains a relevant factor up to inter- 
QD distances of the order of the Kondo screening length 
= XpD/5). The RKKY interaction, instead, decays 
at a typical distance of the order of the Fermi length Af- 
As a function of Cq Eq.© becomes 



5 {r) = D exp 



-1 



[2±C Q (5 ,r)} J n ' 



(10) 



that iteratively converges to the odd Doublet correlation 
energy gain. The minus sign in front of Cq in Eq. i|10|) 
corresponds to the even Doublet energy S e (r), the even 
Doublet is the Kondo structure that develops from \S a ), 
see Eas.l(5|l. H5b|l . Note that Cq couples to the "connec- 
tivity" factor, the 2 in the exponent of Sk, small changes 
in this factor produce exponential changes in S . 
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FIG. 2: Correlation energy gains 5 e and 8 as a function of 
r, for J n = 0.07238 {S K = 0.001D). These energy gains are 
greater than 28 k up to r ~ £x- 

In Fig|3 we plot S e and 8 as function of r. The pe- 
riodicity with which S and S e alternate (~ Xp) is de- 
termined by Cq, S is greater than S e for Cq{t) > 
and viceversa. They cross each other at the tk points 
for which Cq(5k,tk) — 0, and at these points both are 
equal to Sk- For most of the (J„,r) range of interest 
these energies are greater than the RKKY energy 
In turn, these Doublets are connected through the Fermi 
Sea state to form a composite Kondo Singlet with a fur- 
ther correlation energy gain 7 < Sk which is also plotted 
in Fig|5J This second stage Kondo energy has its maxima 
at the zeros Cq, at which 7 = Sk, given for these points 
a total Kondo energy gain (S + 7) equal to 28k- Besides 
these regions 7 is exponentially small. This composite 
singlet is mainly formed using the lower energy Doublet, 
which thus determines the QD-spin correlation. In Table 
[|]we list the main characteristics of the Kondo structures 
present in the two-QD Anderson-Kondo system. 

For these Doublets the QD spin-spin correlation 
(Sl-Sr) is determined by 

(D\S L .S R \D) _ 3 D Q (S,r) 

(D\D) 4 (2±D Q (5,r)) ' 1 > 

where Dq — (OsIq)/ \8sIk) and the upper (lower) sign 
holds for the odd (even) Doublet. 

In Fig[3]it can be seen that the Doublets favor a fer- 
romagnetic alignment of the QD spins. This is due to 
the fact that the enhanced channel is the aa one, a 
ferro configuration fEas. l|4bl5b(l ) and that the 00 chan- 
nel, which is a ferro or antiferro configuration depend- 
ing on the relative phase of their components, is also 



TABLE I: Vertex, bunch states (channels), connectivity fac- 
tors (per channel and total) and total energy for the Kondo 
structures in the two-QD system. From top to bottom: one- 
impurity singlet, odd and even two-impurity Doublets and 
two- impurity composite Super-Singlet. 
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"The channel corresponding to the Doublet of lower energy (higher 
S) dominates the Super-Singlet structure. 




FIG. 3: QD-spin correlation for the doublets. The thicker 
section in each line is the r-region dominated by the cor- 
responding Doublet. The straight segments at 1/4 (—3/4) 
corresponds to the FM (AF) state. 



mainly a ferro configuration at the maxima of S. Clearly, 
the screening action of the hole is more effective if the 
QD-spins are aligned and the distance between them is 
near the resonant condition for the lowest energy holes 
(r ~ riA^/2, |cosfci?r| ~ 1), see Ea . (|5b|l . At the transi- 
tion points from one Doublet to the other it is shown in 
FigE a little negative value for (Sl-Sr) and then a jump 
to positive values. Actually, at those transition points 
the Super Singlet has a maximum in its correlation en- 
ergy (7 ~ Sk) and it is formed with similar weights in 
both Doublets, thus an average of the response of the odd 
and even Doublets is to expect at those regions, given a 
smooth transition with a (Sl-Sr) value near zero from 
one Doublet region to the next one. Note that this mainly 
ferromagnetic QD-QD response is obtained without in- 
clude any possible effect of the RKKY interaction in the 
Doublets structure. 

In FigEl it is shown also the RKKY prediction, based 
upon E R (r) ^ (ErM is equal to 4 In 2 3 2 n (1 - 
2Si(2fcp<r)/7r) D in 1D|2J|). Taken into account both in- 
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teractions, the stronger of them determines the (Sl-Sr) 
response of the system, for max(<5 ,(5 e ) ^> the re- 
sponse is determined by the Doublets curve (Eq.JTTJ, 
Fig|3J) whereas that in the opposite situation the QD- 
spin correlation tends to that of the RKKY. 
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FIG. 4: (Sl-Sr) ID "Quantum Phase Diagram" for the two- 
QD spin correlation. The horizontal dot lines mark the values 
of 5k (in units of D ) for the corresponding values of J n . 

In FigQ] we show the (Sl-Sr) "Parameters Phase 
Diagram" for the ID two-QD system, determined by 
max(<5 ,<5 e ) = |S.r|. This is not a true Quantum Phase 
diagram given that not the Doublets nor the FM (AF) 
states are the ground state of the system, it indicates the 
tendency of the QD-spin correlation to follow a Doublet- 
like or RKKY-like behavior. Over the near vertical 
dashed lines that mark the limits between the odd and 
even Doublets, determined by Cq(8k) = 0, there is lit- 
tle enhancement of the Doublets energies. These zones 
coincide with the AF maxima of the RKKY interaction 
and thus over these lines there are the highest indents of 
the RKKY-AF states, which favor an antiferromagnetic 
alignment of the QD spins. Note also that in the RKKY- 
FM regions (i.e. Y.r 3> <5>) the effect of the correspond- 
ing Doublet is to screen the total Spin of the svstem[l3j 
but, simultaneously, it reinforces the ferromagnetic QD- 
spin response. 

In this work we have shown how a first-order corre- 
lation present in the Two-Impurity Anderson Hamilto- 
nian leads to the formation of a pair of Kondo Dou- 
blet states. The localized spins correlation (Sl-Sr) re- 
sponse for these Doublets strongly differs from the one 
predicted by the second-order RKKY interaction (Fig[3J). 
These Kondo Doublets generates a mainly ferromagnetic 
QD-spin correlation without resort to the second-order 
RKKY interaction. The correlation energy gain of these 
Doublets, which is driven by direct coherence effects of 
the hybridization, exceeds that of the RKKY interac- 
tion in most of the Hamiltonian parameter-space ((Jn, r), 



FigQJ . A remarkable property of this interaction is the 
long distances (~ £k) at which the impurities remain 
in a correlated state. This is due to the fact that just 
the electrons whose energy differs from that of the Fermi 
level less than S contribute significatively to this corre- 
lation. Angular decoherence effects reduce this range for 
a higher-Dimension electron bath gas. These properties 
put these Doublets states well in the experimentally ac- 
cessible range for QD systems built on semiconductor 
devices j^j, were the most relevant parameters can be 
controlled by gate voltages. The VWF method we use 
can be extended to analyze systems with several QD. 
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